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Abstract 



The oscillation frequency Amj of B°B° mixing is measured using the partially- 
reconstructed semileptonic decay B° — > £~VD* + X. The data sample was 
collected with the CDF detector at the Fermilab Tevatron collider during 
1992-1995 by triggering on the existence of two lepton candidates in an event, 
and corresponds to about 110 pb" 1 of pp collisions at ^fs = 1.8 TeV. We 
estimate the proper decay time of the B° meson from the measured decay 
length and reconstructed momentum of the £~D* + system. The charge of 
the lepton in the final state identifies the flavor of the B° meson at its decay. 
The second lepton in the event is used to infer the flavor of the B° meson 
at production. We measure the oscillation frequency to be Arrid = 0.516 ± 
0.099 ^0:035 P s_1 ' wnere the first uncertainty is statistical and the second is 
systematic. 



PACS numbers: 14.40.Nd, 13.20.He 
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I. INTRODUCTION 



Particle-antiparticle mixing in the B°B° system has been known for a decade now [|TJ. 
The phenomenon can be understood as a second-order weak interaction effect. The frequency 
of the oscillation between the two states corresponds to the mass difference Am^ between 
the two mass eigenstates of the B°B° system, B^ and B^. It can be calculated @] from box 
diagrams, where contributions of the top quark in the loop are dominant. Measurements of 
B°B° mixing can therefore determine the magnitude of the Kobayashi-Maskawa matrix || 
element V t d- 

Experiments at the T(4S) resonance have measured the probability of mixing, Xd, inte- 
grated over decay time ||. Experiments at LEP HQ and the Collider Detector at Fermilab 
(CDF) [0, where B hadrons are produced at higher energies, examine the time development 
of mixing and measure the oscillation frequency Am^. By now the Am^ measurements, as 
well as the top quark mass measurements, have become sufficiently precise that other uncer- 
tainties, in particular the B° meson decay constant, limit the precision of the extraction of 
\Vu\- 

The same phenomenon of particle-antiparticle oscillations is expected for the B° S B° S sys- 
tem, where the relevant element of the KM matrix is \V ts \- Due to the difference in the 
involved matrix elements, B° S B° S mixing is expected to proceed with a higher oscillation fre- 
quency, and so far only lower limits on the frequency Am s have been placed |6|,§]. Once the 
B° S B Q S oscillation is established, a measurement of the ratio of the two oscillation frequen- 
cies, Am s /Amd, would provide a useful constraint on the ratio of the KM matrix elements 
|Ks|/|Kd| with less theoretical uncertainty. 

In this paper we report a measurement of B°B° mixing using partially reconstructed 
semileptonic decays. The data used in this analysis were collected in 1992-1995 with the 
CDF detector at the Fermilab Tevatron proton-antiproton collider at a center-of-mass energy 
of y/s =1.8 TeV, and correspond to an integrated luminosity of about 110 pb -1 . We use 
a data sample where events are collected on the existence of two lepton candidates. In 
order to identify semileptonic decays of B mesons, we select events with a lepton (e - or 
H~ , denoted by t~) associated with a D* + meson. (Throughout this paper a reference to a 
particular charge state also implies its charge conjugate.) The £~D* + pairs consist mostly of 
B° decays. The D* + decays are reconstructed using the decay mode D* + — > D tt + , followed 
by D° — > K~tt + , K~7t + 7t + 7t~ , or f^~7r + 7r°. About 500 such decays are reconstructed in the 
data sample. We reconstruct their decay vertices and estimate the proper decay length of the 
B° meson using the momentum of the £~D* + system. The charge of the final state lepton 
identifies the flavor of the B° meson at the time of its decay (£~D* + for B°, and £ + D*~ for 
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B°). The B° meson flavor at its production is inferred from the charge of the second lepton 
in the event (b — > B — > i + vX\ assuming that b and b quarks are produced in pairs. Thus, 
in ideal cases, an opposite-sign lepton pair identifies an unmixed decay, and a same-sign 
pair identifies a mixed decay. We examine decay length distributions of opposite-sign and 
same-sign events and extract the oscillation frequency Am^. 



II. CDF DETECTOR AND TRIGGER 

The CDF detector is described in detail elsewhere |§. We describe here only the detector 
components most relevant to this analysis. Inside the 1.4 T solenoid the silicon vertex detec- 
tor (SVX) [|l(J and the central tracking chamber (CTC) provide the tracking and momentum 
analysis of charged particles. The CTC is a cylindrical drift chamber containing 84 mea- 
surement layers. It covers the pseudorapidity interval \rj\ < 1.1, where rj = — ln[tan(#/2)]. 
In CDF, ip is the azimuthal angle, 9 is the polar angle measured from the proton direction, 
and r is the radius from the beam axis (z-axis). The SVX consists of four layers of silicon 
micro-strip detectors located at radii between 2.9 and 7.9 cm from the beam line and provides 
spatial measurements in the r-ip plane with a resolution of 13 fim. It gives a track impact 



parameter resolution of about (13 + AO/pr) /^rn fll^ l, where pr is the transverse momentum 
of the track with respect to the beam axis and measured in GeV/c. The silicon detectors 
extend to ± 25 cm along the z axis, where z is parallel to the proton beam axis. Since 
the vertex distribution for pp collisions has an rms width of ±30 cm along the z direction, 
a substantial fraction of the interactions occurs outside of the SVX coverage; as a result, 
the average geometric acceptance of the SVX is about 60%. The transverse profile of the 
Tevatron beam is circular and has an rms spread along both x and y axes of ~ 35 /im for 
the data taking period in 1992-1993 and ~ 25 /mi in 1994-1995. The px resolution of the 
CTC combined with the SVX is a{p T )/p T = [(0.0066) 2 + (0.0009 p T ) 2 } 1/2 - Electromagnetic 
(CEM) and hadronic (CHA) calorimeters with projective tower geometry are located outside 
the solenoid and cover the pseudorapidity region \r]\ < 1.1, with a segmentation of Aip =15° 
and At] ~ 0.11. A layer of proportional chambers (CES) is embedded near shower maximum 
in the CEM and provides a more precise measurement of electromagnetic shower profiles and 
an additional measurement of pulse height. A layer of proportional chambers (CPR) is also 
installed between the solenoid and the CEM and samples the electromagnetic showers at 
about one radiation length. Two muon subsystems in the central rapidity region (\r]\ < 0.6) 
are used for muon identification: the central muon chambers (CMU) located just behind 
the CHA calorimeter, and the central upgrade muon chambers (CMP) which lie behind an 
additional 60 cm of steel. The central muon extension chambers (CMX), covering a rapidity 
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region up to < 1.0), are also used. 

CDF uses a three-level trigger system, where at the first two levels decisions are made 
with dedicated hardware. The information available at this stage includes energy deposits 
in the CEM and CHA calorimeters, high p T tracks found in CTC by a track processor, and 
track segments found in the muon subsystems. At the third level of the trigger, the event 
selection is based on a version of off-line reconstruction programs optimized for speed. The 
lepton selection criteria used in level 3 are similar to those described in the next Section. 

Events containing semileptonic B decays and used for this analysis are collected using 
two triggers that require two lepton candidates in an event. The first trigger requires both 
an electron candidate and a muon candidate. The E T threshold for the electron is 5 GeV, 
where E T = Esin9, and E is the energy measured in the CEM. In addition, a track is 
required in the CTC with p T > 4.7 GeV/c that points at the calorimeter tower in (p. The 
muon candidate requires a track in the CTC with matched track segments in the CMU or 
CMX system corresponding to a particle with p T > 2.7 GeV/c. The second trigger requires 
two muon candidates, where the pt threshold is 2.2 GeV/c for each muon, and at least one 
of the muons is required to have track segments in both the CMU and CMP chambers. 

III. RECONSTRUCTION OF SEMILEPTONIC DECAYS OF B MESONS 

The analysis starts with identification of lepton candidates. We require at least two good 
lepton candidates in an event. We then look for the charm meson D* + near each lepton 
candidate to identify the B meson decay B — > £~i>D* + X. A proper correlation between the 
lepton charge and the charm flavor, namely i~ with D* + , and not i + with D* + , is required. 
This decay is used to measure the proper decay length of the B° meson and to identify the 
decay flavor. The charge of the other lepton candidate in the event is used to infer the flavor 
of the B° meson at its production. 

A. Lepton identification 

The identification of electrons makes use of information from both calorimeters and track- 
ing chambers. We require the following: 

• Longitudinal shower profile consistent with electrons, i.e., small leakage of energy into 
the CHA. 
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Lateral shower profiles measured with the CEM [|TTJ and the CES (12| consistent with 
test beam data. 



• Association of a high track with the calorimeter shower based on position matching 
and energy-to-momentum ratio. 

• Pulse heights in the CES and CPR consistent with an electron. 

Photon conversion electrons, as well as the Dalitz decays of 7r° mesons, are removed by look- 
ing for oppositely charged tracks that have small opening angles with the electron candidate. 

Muons are identified based on the geometrical match between the track segments in the 
muon chambers and an extrapolated CTC track. We compute the x 2 of the matching, where 
the uncertainty is dominated by multiple Coulomb scattering in the detector material. We 
require \ 2 < 9 in the r-tp view (CMU and CMP) and x 2 < 12 in the r-z view (CMU). For 
muon candidates in the CMX we require \ 2 < 9 in both the r-tp and r-z views. 

B. Charm meson reconstruction 

To identify B — > £~uD* + X candidates, we search for D* + — > D tt + decays in the vicinity 
of a lepton candidate using two fully reconstructed D° decay modes, D° — > K~n + and 
D° — > K~7t + 7t + 7t~ , and one partially reconstructed mode, D° — > K~7t + tt°. To reconstruct 
D° — > K~7i + decays, we first select oppositely charged pairs of particles using CTC tracks, 
where the kaon mass is assigned to the particle with the same charge as the lepton, as 
is the case in semileptonic B decays. The kaon (pion) candidate is then required to have 
transverse momentum above 1.2 (0.4) GeV/c, and to be within a cone of radius AR = 0.8 
(1.0) around the lepton in i]-ip space, where AR = [(At]) 2 + (Ay?) 2 ] 1 / 2 . To ensure accurate 
decay length measurement, each candidate track, as well as the lepton track, is required to 
be reconstructed in the SVX with hits in at least two layers out of the possible four, and 
with x 2 < 6 per hit. To reduce combinatorial background, we require the decay vertex of the 
D° candidate to be positively displaced along its flight direction in the transverse plane with 
respect to the position of the primary vertex. The primary vertex is approximated by the 
position of the Tevatron beam, which has been determined using independent events ||14|| . 
For the D° — > K~7t + 7t + tt~ mode, the kaon (pion) candidate is required to have transverse 
momentum above 1.2 (0.5) GeV/c, and to be within a cone of radius AR = 0.6 (1.0) around 
the lepton candidate. For the D° — > K^tt+tt mode, the kaon (pion) candidate is required 
to have transverse momentum above 1.2 (0.4) GeV/c, and to be within a cone of radius 
AR = 0.7 (0.8) around the lepton candidate. 
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In order to qualify as a candidate for the signal, the D° candidate has to be in the mass 
range 1.83 to 1.90 GeV/c 2 for the fully reconstructed mode D° — ► K~n + and in the range 
1.84 to 1.88 GeV/c 2 for the — > K~tt + 7t + tt~ mode. For the partially reconstructed mode 
D° — > K~7t + 7t , we require the mass of a K~tt + pair to be between 1.5 and 1.7 GeV/c 2 ; 
we do not reconstruct the 7r° meson and in the subsequent analysis treat the K~tt + pair 
as if it were a D° meson. For each mode, we reconstruct the D* + meson by combining an 
additional CTC track, assumed to have the pion mass, with the D° candidate, and computing 
the mass difference, AM, between the D°7r + and D° candidates. Figure p] shows the AM 
distributions for the three D° decay modes. In Fig. [l|(c) the peak is broadened because 
of the missing 7r° meson. The dotted histograms show the spectra from the "wrong sign" 
(D°7T~) combinations, where no significant signals are observed. We define the signal region 
as follows: the two fully reconstructed modes use the AM range 0.144 to 0.147 GeV/c 2 , 
and the K~tt + it mode uses the range AM < 0.155 GeV/c 2 . The numbers of events in 
the signal regions are 216, 256, and 416 for the three modes. We estimate the numbers of 
combinatorial background events by using the shapes of the AM spectra of the wrong sign 
(D°7r~) combinations and normalizing them to the number of events in the AM sideband. 
The estimated background fractions are 0.227 ± 0.036, 0.326 ± 0.040 and 0.543 ± 0.050, 
respectively. They are summarized in Table |. 



C. Sample composition 

Apart from combinatorial backgrounds, the £~D* + signal sample contains events which 
originated from physics sources other than the B° meson decays. The main contribution 
comes from B~ meson decays. The semileptonic decays of B mesons can be expressed as 
B — > £~VD, where D is a charm system whose charge is correlated with the B meson 
charge. If only the two lowest mass charm states, pseudoscalar (D) and vector (D*) mesons, 
are produced, the £~ D* + combination can arise only from the B° decay. However, it has been 
known that the above two lowest mass states do not saturate the total semileptonic decay 
rates. All data indicate that higher mass charm mesons, D** states, as well as non-resonant 
D^tx pairs, are responsible for the rest of the semileptonic decays [13[]. In this analysis we 



do not distinguish resonant and non-resonant components, and refer to both of them as D** 
mesons. 

These D** meson decays can dilute the charge correlation between the final states and 
the parent B meson. For example, the D**° meson can be produced by the decay B~ — > 
£~vD**°, which subsequently can produce both D* + n~ and D* n final states. This results 
in misidentification of the B~ meson decay as B° -> D* + £-vX. Nevertheless, the £~D*+ 
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combination is dominated by -B meson decays. 

In order to estimate the fraction g~ of B~ decays relative to the sum of B~ and B° 
mesons in the observed £~D* + sample, we follow the method used in the CDF measurement 
of the B~ and B° meson lifetimes jyj using semileptonic decays. We describe the method 
here as well. 

The production rates of charged and neutral B mesons and their semileptonic decay 
widths are assumed to be equal. We also assume the D** mesons decay exclusively to a D^ir 
pair via the strong interaction, thereby allowing us to determine the branching fractions, 
e.g. D(*> + 7r vs. D^°tc + , using isospin symmetry. We consider three factors affecting the 
composition. First, the composition depends on the fraction /** of the D** mesons produced 
in semileptonic B decays, 

= B(B^£-uD**) _ 1 B(B -> t-vD) + B(B -> i~vD*) 



B(B-^e-uDX) B(B^i-uDX) 

where B denotes a branching fraction and B is a B~ or B° meson. The CLEO experiment 
measures the fraction of exclusive decays to the two lowest mass states to be 0.64 ± 0.10 ± 



0.06 [[[(J. Thus, we estimate that /** = 0.36 ± 0.12. A few experiments have recently 
observed some D** modes [[H]], but the sum of exclusive modes still does not equal the total 
semileptonic rate. Second, the fraction g~ depends on the relative abundance of various 
possible D** states, because some of them decay only to D*ir and others to Dtt, depending 
on the spin and parity. The abundance is not measured very well at present. Changing the 
abundance is equivalent to changing the branching fractions for D*7r and modes averaged 
over various D** states. We define the quantity 

_ B(D** -> D*n) 

v = B(D** -f D*n) + B(D** -f Dti) ' 

We assume the relative abundance of the four D** mesons predicted by the Isgur-Scora- 
Grinstein-Wise (ISGW) model [|18|], which corresponds to Py = 0.64. After inclusion of 
non-resonant contributions, we use Py = 0.65 as our nominal choice. We also consider the 
values Py = 0.26 and 1.00. Third, the composition depends on the ratio of the B~ and 
B meson lifetimes, because the number of £~D* + events is proportional to the semileptonic 
branching fraction, which is the product of the lifetime and the partial width. We use the 
ratio t(B-)/t(B°) = 1.02 ± 0.05 . 



We also take into account the differences in the reconstruction efficiencies for the B — > 
£~i?D* and D** decay modes. We examine this effect by using Monte Carlo events where 
the ISGW model is used to describe the semileptonic decays. We shall describe the Monte 
Carlo simulation later. We find that the D** mode has an efficiency that is lower than in 
the D* mode by about 50% (25%) for leptons above 5 GeV/c (2 GeV/c). 
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We find that g~ = 0.19±g;?g for the yr D* + sample and g~ = 0.14 ±g;g| for the e~D*+ 
sample. The central values correspond to the nominal choice of the parameters, /** = 0.36, 
Py = 0.65, and t(B~)/t(B°) = 1.02. The uncertainties reflect maximum changes in g~ 
when /**, Py and the lifetime ratio are changed within their uncertainties, namely /** to 
0.24 and 0.48, Py to 0.26 and 1.0, and t(B-)/t(B°) to 0.97 and 1.07. The difference between 
the muon and electron channels arises from the difference in kinematic requirements. 

There are other physics processes that can produce the lepton-D* + signature. The 
largest background comes from the decay of the B° s meson, B° s — > £~T7D** + , followed by 
D** + — > D* + K°. This process is estimated to contribute about 3% of the lepton-/}* + signal. 
Other processes such as B — > t~v t D* + X followed by r~ — > t~V(V T ^ and B — > D~D* + X fol- 
lowed by D~ — > £~X, are suppressed severely because of branching fractions and kinematic 
requirements on leptons. We ignore these backgrounds here. Therefore, the fraction of B° 
mesons is given by g° = 1 — g~. We treat effects of the physics backgrounds as a systematic 
uncertainty. 



IV. DECAY LENGTH MEASUREMENT AND MOMENTUM ESTIMATE 

The B meson decay vertex Vb is obtained by intersecting the trajectory of the lepton 
track with the flight path of the D° candidate. The B decay length Lb is defined as the 
displacement of Vb from the primary vertex Vp, measured in the plane perpendicular to the 
beam axis, and projected onto the transverse momentum vector of the lepton-D* + system: 



T _ (V B -Vp)-p£ 



Pt 

A schematic representation of the B° meson semileptonic decay is illustrated in Fig. |2|. 

To measure the proper decay length of a B meson, we need to know the momentum of 
the B meson. In semileptonic decays, the B meson momentum cannot be measured precisely 
because of the missing neutrino. We use the transverse momentum of the observed system, 
p l T D * + , to estimate the B meson transverse momentum for each event. We denote the 
ratio of the two momenta by K = p^ D * + /pt, an d introduce a corrected decay length defined 
as 

- T mB IV\ 

x = L B e _ Dt+ {K), 
Pt 

which we call the "pseudo-proper decay length." The average correction for the missing mo- 
mentum is achieved by the constant (K) . The correction for a finite width of the distribution 
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of the ratio K is performed during fits to decay length distributions. We shall describe the 
fits later. 

A typical resolution on this decay length x due to vertex determination is 50 /xm, including 
the contribution from the finite size of the primary vertex. For subsequent decay length 
measurements, we use only those events in which the resolutions on reconstructed decay 
lengths are smaller than 0.05 cm. We also require the proper decay length of the D° meson, 
measured from the B meson decay vertex to the D° decay vertex, to be in the range from 
—0.1 cm to 0.1 cm, with its uncertainty smaller than 0.05 cm. These cuts reject poorly 
measured decays and reduce random track combinations. In addition, we limit ourselves to 
events with reconstructed decay lengths in the range between —0.15 cm and 0.3 cm. These 
cuts have been applied already for the charm signals shown in Fig. [l]. 

The distribution of the momentum ratio K is obtained from a Monte Carlo calculation. 
The b quarks are generated according to the px spectrum by the QCD calculation in the 
next-to- leading order |3U|| . The fragmentation model by Peterson and others [21] is used. 
The CLEO event generator |Q is then used to describe the B meson decays. In particular, 
the semileptonic decays adopt the ISGW model flT8| . A typical K distribution thus obtained 
has an average value of 0.85 with an rms width of 0.14, and shows only a weak dependence 
on p T D * + in the range of interest, which is typically between 10 and 20 GeV/c 2 . It is 
also independent of the D° decay mode except for the partially reconstructed mode D° — > 
K^7i + 7T , which has a slightly lower mean value (about 0.80) because of the missing 7r° 
meson. Typical K distributions are shown in Fig. |3|. 

We fit the observed pseudo-proper decay length distributions for both opposite-sign and 
same-sign events. This fit determines parameters that will be used later in the fit for the 
oscillation frequency Am^. It also yields the B° meson lifetime as a check of the momentum 
correction described above. We use the maximum likelihood method. The likelihood used 
to fit the events in the signal region is expressed as 

£sig = ITK 1 — /bg)-^sig(^) + /bg-^bg(^)], 

i 

where x, is the pseudo-proper decay length measured for event i, and the product is taken 
over observed events in the sample. The first term in the likelihood function represents 
the contribution of B decay signal events, while the second term accounts for combinatorial 
background events whose fraction in the sample is /bg- 

The signal probability density function J-'sidx) has two components and is expressed as 

Fsig(x) = S^sIgO) + (1 - £~)^igO), 

where F^idx) and F^idx) are the normalized probability density functions for the B~ and 
B° meson decays, respectively, and g~ is the fraction of B~ mesons as defined earlier. Each 
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component consists of an exponential decay function, denned for positive decay lengths, 
smeared with a normalized K distribution D{K) and a Gaussian distribution with width 
serf. 

= JdKD->°(K) 

where r is the appropriate B meson lifetime, c is the speed of light, 6{x) is the step function 
defined as 6(x) = 1 for x > and 9{x) = for x < 0, and the symbol "(g)" denotes a 
convolution. G(x) is the Gaussian distribution given by 

G(x) = i exp ("i) • 

where <jj is the estimated resolution on Xj. The scale factor s is introduced in order to 
account for a possible incompleteness of our estimate of the decay length resolution. The 
integration over the momentum ratio K is approximated by a finite sum 

J dKD(K) -> ^2 D(Kj)AK, 

j 

where the sum is taken over bin j of a histogrammed distribution D(Kj) with bin width 
AK = 0.02. The K distributions for B~ and B° mesons are slightly different because the 
B~ — > £~i?D* + X decay involves more missing particles. 

The pseudo-proper decay length distribution of combinatorial background events, 
jF BG (x), is measured using AM sideband events, assuming that they represent the com- 
binatorial background events under the signal mass peaks. The functional form of the dis- 
tribution is parameterized empirically by a sum of a Gaussian distribution centered at zero, 
and positive and negative exponential tails smeared with a Gaussian distribution: 

F BG (x) = (1 - /_ - /+) G(x) 

+ J-e^ x ) exp ^ + ^j ®G(x). (2) 

The shape of the background function (parameters f± and \±) and the resolution scale 
factor s, as well as the B meson lifetime cr, are determined from a simultaneous fit to signal 
and sideband events. We fix the ratio of the B~ and B° meson lifetimes and fit for the 
B° meson lifetime only. To determine those parameters, we use the combined likelihood 
C defined as C = £sig £bg, wriere £bg = Ilfc^BG^fc) an d the product is taken over 
event k in the background sample. The amount of combinatorial background /bg is also a 



n, n K ( Kx \ „t s 

9{X) ^)^[-^K))® G{X) ' 



14 



parameter in the simultaneous fit. This parameter is constrained by adding a term = 
|(/bg — (/bg)) 2 /°"bg ^° the negative log-likelihood — £ = — ln£. The average background 
fraction (/bg) an d its uncertainty <tbg axe estimated from the signal mass distributions and 
are given in Table |. 

The background sample for the £~D* + candidates is taken from the AM sidebands: we 
use the right sign (D°7r + ) sideband 0.15 < AM < 0.19 GeV/c 2 for the two fully reconstructed 
D° modes, and 0.16 < AM < 0.19 GeV/c 2 for the D° -> R-n+n mode. We also use the 
wrong sign pion combinations in the range AM < 0.19 GeV/c 2 for all three D° decay modes. 
The background samples are summarized in Table [nj. 

The pseudo-proper decay length distributions of the background samples are shown in 
Fig. |], together with fit results. The background parameter values and the resolution scale s 
determined from the fit are listed in Table |TTT| . The corresponding decay length distributions 
of the signal samples are shown in Fig. We find the lifetimes to be ct(B°) = 470 ± 
44, 407 ± 40 and 419 ± 39 [xm for the three D° decay modes. The quoted uncertainties are 
statistical only. When a combined fit to the three modes is made, we find ct(B°) = 433 ±24 
/im. These results are consistent with the world average value of 468 ± 12 ^m [13|| . 



V. B°B° MIXING MEASUREMENT 



The probability that a B° meson at t = decays as -B (unmixed) or as B° (mixed) at 
a proper time t is given by 

-Ptjnm (t) = 7^I ex P (~~) (! + cos Am d t), 

Puix(t) = ^I ex P (~~) ^ ~ cosAm d t), 

where r is the B° meson lifetime, and we have ignored CP violation and the width difference 
Ar between the two mass eigenstates of the B°B° system. We determine the mixing pa- 
rameter Amd by a simultaneous fit to the decay length distributions of unmixed and mixed 
decay events. 

We have reconstructed the B° meson decay B° — > £~i?D* + X . The charge of the lepton 
identifies the flavor of the B° meson at its decay. In order to infer the B° meson flavor at 
its production, we use the second lepton candidate, which is presumed to originate from the 
other B hadron in the event. When the other B hadron, containing the b quark, decays 
semileptonically, it produces a positively charged lepton £ + . If the B° meson, reconstructed 
in the £^uD* + X decay mode, decayed in an unmixed state, the two leptons in the event 
would have the opposite charge. Similarly, if the B° meson decayed in a mixed state, the 
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leptons would have the same charge. Therefore, in the ideal case, the opposite-sign (OS) 
events identify the unmixed decays of the B° meson, while the same-sign (SS) events identify 
the mixed decays. However, the second lepton can originate from the sequential decay of B 
hadrons, b — » c — > £~X, or from a mixed decay of the neutral B mesons, h — > B° (B®) — > B° 
(B®) — > £~X. The lepton candidate could also be a misidentified hadron. In these cases the 
second lepton candidate will not identify the production flavor correctly. In order to account 
for these possibilities, we introduce the probability of flavor misidentification and denote it 
by W. 

As mentioned above, we classify events depending on the sign (OS or SS) of the two lepton 
candidates in an event. A finite flavor misidentification probability W results in moving 
unmixed decay to the same-sign sample and mixed decays to the opposite-sign sample. 
Thus, we obtain the following probability distributions for the opposite-sign and same-sign 
events: 

P os (t) = (1 - W0P UNM (t) + WP Mlx (t) = i- exp (-^j [1 + (1 - 2W) cos Am d t], 
P ss (t) = (1 - W)P Mlx (t) + WP v Mt) = y t exp (-^) [1 - (1 - 2W) cos Am d t}. 

From these expressions it is evident that the flavor misidentification probability does not 
affect the oscillation frequency, although it does reduce its amplitude by a factor 1 — 2W. 
We determine the two quantities that appear in the above expression, Am^ and W, simul- 
taneously from the data sample by examining the decay length distributions. 



A. Am c i fit 

We use the maximum likelihood method to extract the oscillation frequency Am^. The 
likelihood is given by C = Ili where Xi is the pseudo-proper decay length measured 

for event i, and the product is taken over events in the signal sample. The likelihood function 
T(x) is expressed as follows, depending on the sign (OS or SS) of an event: 

nx) = [ (1 - /bg) ^ G {x) + / BG (1 - /ss) ^bg(x) if OS, 
"\(1- /bg) ^s S ig(^) + /bg/ss ^bg(x) if SS, 

where /bg is the fraction of combinatorial background events in the sample. The background 
function J-bg(^) is the same as in the lifetime fit (Eq. (|2[)), and its shape is taken to be the 
same for both opposite-sign and same-sign events, /ss is the fraction of same-sign events in 
the combinatorial background. 

Each of the signal functions consists of two components, one for the B° meson and the 
other for the B~ meson: 
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*skU*) = (1 - 9-) [(1-W) ^ NM (x) + W J* lx (x) ]+ g -(l-W) F-{x), 
HUx) = (1 - 9~) [(1-W) Kix(x) + WF° NM (x) ] + g-WT-(x\ 

where g~ is the fraction of B~ meson decays among the signal, and W is the flavor misiden- 
tification probability of the second lepton. The B° component of the opposite-sign function 
jFg^(x) contains two terms: the first term represents correctly tagged unmixed decays (prob- 
ability 1 — W) , while the second term represents incorrectly tagged mixed decays (probability 
W). Similarly, the B° component of the same-sign function J^siqIx) consists of correctly 
tagged mixed decays and incorrectly tagged unmixed decays. Since the B~ meson does not 
mix, it appears in the opposite-sign function when the production flavor is tagged correctly, 
and in the same-sign function when tagged incorrectly. The B~ function J r ~(x) is a smeared 
exponential decay function and is the same as in the lifetime fit (Eq. (0)). The B° functions 
have an additional factor for the mixing and are given by 



?VNM,mx( x ) = J dKD(K) J fl(x) exp 



K Kx 



(Am d K > 

1 ± COS TT7T X 

\ C \ K ) J 



®G{x) 



2ct(K) * \ ct(K}^ 

where the sign + (— ) before the cosine corresponds to the unmixed (mixed) decay function. 
The background and B~ functions are normalized so as to give unity when integrated over 
x. The B° functions give unity when integrated over x and summed over the two decay 
possibilities, unmixed and mixed. The free parameters in the fit are the oscillation frequency 
Amd and the flavor misidentification probability W. We fix the shape of the background 
function (parameters f± and X±) and the resolution scale factor as determined from the 
background sample. The lifetime of the B° meson is also fixed to the value determined 
earlier in the signal sample. This procedure has been found [^J to improve slightly the 
sensitivity in Amd determination. It is confirmed with our study using Monte Carlo events. 
The background fraction /bg m the sample and the same-sign fraction /gg of the background 
are fit parameters, but they are constrained by adding a x 2 term to the negative log-likelihood 
— i = — In C, as in the lifetime fits. 

The B meson pseudo-proper decay length distributions of the £~D* + signal sample are 
shown in Fig. |]for the opposite-sign and same-sign events, as well as for the sum of the two. 
The three D° decay modes are combined there. We find 498 opposite-sign events and 390 
same-sign events. The fraction of same-sign events in the combinatorial background, /gg, is 
estimated using events in the background sample. The estimated same-sign fractions fss are 
summarized in Table [TV]. The fit results are Am d = 0.516±0.099 ps _1 and W = 0.325±0.033, 
where uncertainties are statistical only. They are summarized in Table [V] along with other 
fit parameters. The fit results are also shown in Fig. [| The oscillatory behavior can be 
seen more directly when the asymmetry between the opposite-sign and same-sign events is 
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examined as a function of the pseudo-proper decay length. We define the asymmetry as 

N os (x) - N ss (x) 
^ ~ iV os (x) + iV ss (x)' 

where N os (x) (N ss (x)) is the number of opposite-sign (same-sign) events. In the ideal 
case where the backgrounds, the flavor misidentification, and the decay length smearing are 
absent, the asymmetry is given by A(x) = cos(Am^x/c). The asymmetry distribution of 
the signal sample is illustrated in Fig. [7|, together with the fit result. 

B. Systematic uncertainties 

The sample composition is a source of systematic uncertainty in the oscillation frequency 
measurement. We have described it in terms of the parameters /**, Py and the lifetime 
ratio t(B~)/t(B°). We change each one of the parameters to another value while keeping 
the others at their nominal values, compute the sample composition g~ , and repeat the fit 
procedure for Am^. We note that the momentum correction factors (K distributions) need 
to be modified accordingly; the K distributions for the decay B — > l~i)D** have lower mean 
values because of additional missing particle(s), and changing the amount of D** decays 
results in changes in the K distributions. The results are summarized in Table [VT. We 
interpret the observed changes as systematic uncertainties. 

Other sources of systematic uncertainties considered in this analysis are summarized 
in Table |V11| . Physics background processes are studied by adding their simulated decay 
length distributions to the background function. In addition, the shapes of the decay length 
distributions of the combinatorial background events and the signal lifetime are subject to 
uncertainty because they are determined with finite statistical precision. They are changed 
within uncertainties, and the fit is repeated. We interpret the observed changes as the 
systematic uncertainty due to this source. 

Other sources of systematic uncertainties include our estimates of the decay length res- 
olution and of the B meson momentum. We have introduced a resolution scale factor s and 
find a value of about 1.2. We change this factor to 1.0 or 1.4 and repeat the fit. We assign 
the observed changes as an uncertainty. The B° meson momentum estimate (K distribution) 
is subject to some uncertainty too, because it depends on the kinematics of B meson pro- 
duction and of semileptonic decays. We investigate different production and decay models 
using the procedure described in Ref. ||15|| , and estimate the uncertainty in the B° meson 



momentum to be 2%, which translates directly to the Ama uncertainty. 

All contributions are added in quadrature to give the total systematic uncertainty in 
Arrid of taois P s_1 an< ^ ^ n ^ ne fl avor misidentification probability W of I o!o?2- 
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VI. CONCLUSION 



We have measured B°B° mixing using the semileptonic decay B° — > l~i)D* + X recon- 
structed among pp collision events with two lepton candidates. The proper decay length is 
estimated from reconstructed decay vertices and the momentum of the £~D* + system. A 
high B° purity and a relatively good momentum resolution are achieved. The second lepton 
candidate in the event is used to infer the flavor of the B° meson at the time of its produc- 
tion, with a flavor misidentification probability of W = 0.325 ± 0.033 io!o?2- The frequency 
of the oscillation is measured to be 

Am d = 0.516 ± 0.099 ±g;g|§ ps~\ 

where the first uncertainty is statistical and the second is systematic. The result is consistent 
with other recent measurements f^]. 

The method could be also applied in the future to a search for B®B° S oscillations with a 
modest value of Am s by reconstructing the meson produced in the semileptonic decay 
B° s -> £-pD+X. 



ACKNOWLEDGMENTS 



We thank the Fermilab staff and the technical staffs of the participating institutions for 
their vital contributions. This work was supported by the U.S. Department of Energy and 
National Science Foundation; the Italian Istituto Nazionale di Fisica Nucleare; the Ministry 
of Education, Science and Culture of Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the Republic of China; the A. P. Sloan 
Foundation; and the Alexander von Humboldt-Stiftung. 



19 



REFERENCES 

[1] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 192, 245 (1987); CLEO Col- 
laboration, M. Artuso et al, Phys. Rev. Lett. 62, 2233 (1989). 

[2] T. Inami and C. S. Lim, Prog. Theor. Phys. 65, 297 (1981); 65, 1772(E) (1981). 

[3] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[4] ARGUS Collaboration, H. Albrecht et al, Z. Phys. C 55, 357 (1992); CLEO Collab- 
oration, J. Bartelt et al, Phys. Rev. Lett. 71, 1680 (1993); ARGUS Collaboration, 
H. Albrecht et al, Phys. Lett. B 324, 249 (1994). 

[5] L3 Collaboration, M. Acciarri et al, Eur. Phys. J. C 5, 195 (1998); DELPHI Collabo- 
ration, P. Abreu et al, Z. Phys. C 76, 579 (1997); OPAL Collaboration, K. Ackerstaff 
et al, Z. Phys. C 76, 417 (1997); ALEPH Collaboration, D. Buskulic et al, Z. Phys. C 
75, 397 (1997); OPAL Collaboration, G. Alexander et al, Z. Phys. C 72, 377 (1996). 

[6] OPAL Collaboration, K. Ackerstaff et al, Z. Phys. C 76, 401 (1997). 

[7] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 80, 2057 (1998); CDF Collaboration, 
F. Abe et al, Phys. Rev. D 59, 032001 (1999); CDF Collaboration, F. Abe et al, Phys. 
Rev. D 60, 051101 (1999); CDF Collaboration, F. Abe et al, FERMILAB-Pub-99/019-E 
(hep-ex/990311) (1999), to be publihed in Phys. Rev. D. 

[8] ALEPH Collaboration, R. Barate et al, Eur. Phys. J. C 4, 367 (1998); DELPHI Col- 
laboration, W. Adam et al, Phys. Lett. B 414, 382 (1997); CDF Collaboration, F. Abe 
et al, Phys. Rev. Lett. 82, 3576 (1999). 

[9] CDF Collaboration, F. Abe et al, Nucl. Instrum. Methods Phys. Res. A 271, 387 (1988), 
and references therein. 

[10] D. Amidei et al, Nucl. Instrum. Methods Phys. Res. A 350, 73 (1994); P. Azzi et al, 

Nucl. Instrum. Methods Phys. Res. A 360, 137 (1995). 
[11] CDF Collaboration, F. Abe et al, Phys. Rev. D 52, 2624 (1995). 

[12] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 68, 2734 (1992); CDF Collaboration, 

F. Abe et al, Phys. Rev. D 48, 2998 (1993). 
[13] C. Caso et al, Eur. Phys. J. C 3, 1 (1998). 

[14] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 71, 3421 (1993); CDF Collaboration, 
F. Abe et al, Phys. Rev. D 57, 5382 (1998). 

[15] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 76, 4462 (1996); CDF Collaboration, 
F. Abe et al, Phys. Rev. D 58, 092002 (1998). 

[16] CLEO Collaboration, R. Fulton et al, Phys. Rev. D 43, 651 (1991). 

[17] ALEPH Collaboration, D. Busklic et al, Phys. Lett. B 345, 103 (1995); DELPHI Col- 
laboration, P. Abreu et al, Z. Phys. C 71, 539 (1996); OPAL Collaboration, R. Akers 
et al, Z. Phys. C 67, 57 (1995); CLEO Collaboration, Phys. Rev. Lett. 80, 4127 (1998). 



20 



[18] N. Isgur, D. Scora, B. Grinstein, and M. Wise, Phys. Rev. D 39, 799 (1989); D. Scora 

and N. Isgur, Phys. Rev. D 52, 2783 (1995). 
[19] R. M. Barnett et ai, Phys. Rev. D 54, 1 (1996). 

[20] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B327, 49 (1989); B335, 260(E) 
(1990). 

[21] C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, Phys. Rev. D 27, 105 (1983). 
[22] P. Avery, K. Read, and G. Trahern, Report No. CSN-212, 1985 (unpublished). 
[23] ALEPH Collaboration, D. Buskulic et ai, Phys. Lett. B 377, 205 (1996). 



21 



TABLES 

TABLE I. Definition of signal samples, numbers of candidates and estimated background 
fractions. 



B 


mode 


D° 


mode 


D° mass range 


AM range 


Events 


Background fraction 










(GeV/c 2 ) 


(GeV/c 2 ) 






£- 


D* + 


K~ 


■7T+ 


1.83 - 1.90 


0.144 - 0.147 


216 


0.227 ±0.036 


t 


D*+ 


K~ 




1.84 - 1.88 


0.144 - 0.147 


256 


0.326 ± 0.040 


r 


D*+ 


K~ 


-vr+vr 


1.50 - 1.70 


< 0.155 


416 


0.543 ± 0.050 



TABLE II. Definition of background samples and numbers of events. 



B mode 


D° mode 


D° mass range 


AM ran§ 


re (GeV/c 2 ) 


Events 






(GeV/c 2 ) 




D»TT- 




e~D*+ 


R~TT + 


1.83 - 1.90 


0.15 - 0.19 


< 0.19 


2418 


£~D* + 




1.84 - 1.88 


0.15 - 0.19 


< 0.19 


5139 


£-D* + 


K-TT+TT 


1.50 - 1.70 


0.16 - 0.19 


< 0.19 


1663 



TABLE III. Background shapes obtained from a simultaneous fit to signal and background 
samples. 



B mode D° mode scale s /+ A + (/im) /_ A_ (//m) 

£~D* + K--K+ 1.21 ±0.05 0.361 ±0.015 474 ± 20 0.157 ±0.015 392 ± 49 

t~D* + R-tt+tt+it- 1.17 ±0.03 0.332 ± 0.012 331 ± 11 0.098 ± 0.010 230 ± 22 

i~D* + R-tt+tt 1.21 ±0.04 0.367 ±0.016 433 ± 22 0.095 ± 0.014 293 ± 33 



TABLE IV. Numbers of opposite-sign (OS) and same-sign (SS) events in the signal and back- 
ground samples, and the fraction /ss of same-sign events in the combinatorial background. 



B mode 


D° mode 




Signal 






Background 


Sum 


OS 


SS 


OS 


SS /ss 


e~D*+ 


K--K+ 


216 


121 


95 


1240 


1178 0.487 ±0.010 


£-D*+ 




256 


146 


110 


2501 


2638 0.513 ± 0.007 


£-D* + 


R-TT+TT 


416 


231 


185 


902 


761 0.458 ±0.012 


£'D* + 


Total 


888 


498 


390 
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TABLE V. Results of the Am d fit. 



Parameter Input Output 

Am d 0.516 ±0.099 

W 0.325 ± 0.033 

/ BG (D° -> K-TT+) 0.227 ±0.036 0.218 ± 0.030 

/ BG (D° -► K-tt+tt+tt-) 0.326 ± 0.040 0.355 ± 0.032 

/ BG (-D K-tt+tt ) 0.543 ± 0.050 0.489 ± 0.034 

fss (D° ^ K~tt + ) 0.487 ±0.010 0.488 ± 0.010 

fss {D° -> K-tt+tt+tt-) 0.513 ±0.007 0.513 ± 0.007 

fss (D° K'tt+tt ) 0.458 ± 0.012 0.465 ± 0.011 



TABLE VI. Measurement of Am d under various sample composition conditions. Quoted 
uncertainties are statistical only. 



J ** 


Pv 


r(B-) 


a' 




Am d 


W 






r(B0) 




e~D* + 


(ps- 1 ) 




0.24 


0.65 


1.02 


0.121 


0.087 


0.497 ± 0.093 


0.323 ± 0.033 


0.36 


0.65 


1.02 


0.187 


0.138 


0.516 ±0.099 


0.325 ±0.033 


0.48 


0.65 


1.02 


0.263 


0.202 


0.536 ±0.108 


0.327 ±0.033 


0.36 


0.26 


1.02 


0.090 


0.063 


0.488 ± 0.090 


0.323 ± 0.033 


0.36 


1.00 


1.02 


0.250 


0.190 


0.532 ±0.106 


0.327 ±0.033 


0.36 


0.65 


0.97 


0.179 


0.132 


0.511 ±0.098 


0.324 ±0.033 


0.36 


0.65 


1.07 


0.194 


0.144 


0.520 ±0.101 


0.325 ± 0.033 



TABLE VII. A summary of systematic uncertainties in Amd measurement. 



Source Contribution to 





Am d (ps x ) 


w 


Sample composition 






D** fraction (/**) 


+0.020 
-0.018 


+0.003 
-0.001 


D** composition {Py) 


+0.017 
-0.027 


+0.002 
-0.001 


Lifetime ratio t{B~)/t(B q ) 


± 0.005 


± 0.001 


Physics background 


+0.003 
-0.000 


+0.000 
-0.010 


Background shape, B° lifetime 


± 0.001 


+0.003 
-0.007 


Decay length resolution 


± 0.005 


+0.003 
-0.002 


B meson momentum estimate 


± 0.011 




Total 


+0.029 
-0.035 


+0.006 
-0.012 
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FIG. 1. Reconstructed D* + — > D°7r + signals in events with two lepton candidates. The D* + 
meson is associated with a lepton (£~) candidate. Distributions of AM for three D° decay modes 
are shown: (a) D° — > K _ 7r + , (b) D° — ► i^ _ 7r + 7r + 7r _ , and (c) Z)° — ► i^ _ 7r + 7r°. Dotted histograms 
show the distributions for wrong sign (D°ir~) combinations. 
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FIG. 2. Schematic representation of the decay B° — > £ uD* + , followed by D* + — ► D°7r + and 
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FIG. 3. Distributions of the momentum ratio K (see text) obtained from Monte Carlo calcula- 
tions for decays B° — > £~ vD* + X , followed by D* + — ► D°7r + . Three D° decay modes are shown: 
(a) D° -► ^-vr+, (b) L>° -► K-vr+^+^-, and (c) L>° -► R-Tt+Tt . 
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FIG. 4. Distributions of B meson pseudo-proper decay lengths for I D* + background samples 
(points). Three decay modes are shown: (a) D* + — > Z)°7r + , D° — ► ivT _ 7r + , (b) Z)* + — > D°7r + , 
L>° -► i^-vr+vr+vr-, and (c) £>*+ -► D°vr+, L>° -► R-ir+ir . Curves show fit results. 
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FIG. 5. Distributions of B meson pseudo-proper decay lengths for l~D* + signal samples 
(points). Three decay modes are shown: (a) D* + — > Z)°7r + , D° — ► -ftT~7r + , (b) Z)* + — > D°7r + , 
£)° — ► ivr _ 7r + 7r + 7r _ , and (c) Z)* + — ► Z)°7r + , D° — > .ftr~7r + 7r . The three modes are combined in (d). 
Also shown are the results of lifetime fits: the B° component (dashed curve), the B~ component 
(dot-dashed curve), the background component (dotted curve), and the sum of all components 
(solid curve). 
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Pseudo — Proper Decay Length (cm) 



FIG. 6. B meson pseudo-proper decay length distribution (points) estimated from the £~D* + 
candidates for (a) opposite-sign events, and (b) same-sign events, and (c) the sum of the two. 
Curves show the result of the Am^ fit: the B° component (dashed curve), the B~ component 
(dot-dashed curve), the background component (dotted curve), and the sum of all components 
(solid curve). 
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FIG. 7. Charge asymmetry of the I D* + candidates as a function of pseudo-proper decay 
length (points). The solid curve shows the result of the Am^ fit. 
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